Accepted for publication to ApJ 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



O 

(N 

o 

Q 

00 



X 
6 



> 
o 

00 

O 



MAGNETIC FIELD DECAY DUE TO THE WAVE-PARTICLE RESONANCES IN THE OUTER CRUST OF 

THE NEUTRON STAR 

HiROYUKi R. Takahashi^, Kei Kotake'^'^, and Nobutoshi Yasutake^ 

Accepted for publication to ApJ 

ABSTRACT 

Bearing in mind the application to the outer crust of the neutron stars (NSs), we investigate the 
magnetic field decay by means of the fully relativistic Particle- In-Cell simulations. Numerical compu- 
tations are carried out in 2-dimensions, in which the initial magnetic fields are set to be composed both 
of the uniform magnetic fields that model the global fields penetrating the NS and of the turbulent 
magnetic fields that would be originated from the Hall cascade of the large-scale turbulence. Our re- 
sults show that the whistler cascade of the turbulence transports the magnetic energy preferentially in 
the direction perpendicular to the uniform magnetic fields. It is also found that the distribution func- 
tion of electrons becomes anisotropic because electrons with lower energies are predominantly heated 
in the direction parallel to the uniform magnetic fields due to the Landau resonance, while electrons 
with higher energies are heated mainly by the cyclotron resonance that makes the distribution func- 
tion isotropic for the high energy tails. Furthermore we point out that the degree of anisotropy takes 
maximum as a function of the initial turbulent magnetic energy. As an alternative to the conventional 
ohmic dissipation, we propose that the magnetic fields in the outer crust of NSs, cascading down to 
the electron inertial scale via the whistler turbulence, would decay predominantly by the dissipation 
processes through the Landau damping and the cyclotron resonance. 
Subject headings: stars: neutron - plasmas - turbulence 



1. INTRODUCTION 

Pushed by the accumulating observations of radio pul- 
sars, accreting neutron stars (NSs), and most recently, 
magnetars, extensive studies have been carried out to un- 
derstand the e volution of magnetic fields i n neutron stars 
(e.g.. Bhattac harya fc SrinivascinI Il995t iHarding fc Lail 
I2006t rReiseneggeril2009D . The known radio pulsars so far 
are generally categorized into two classes: young (< 10'' 
yr) pulsars with spin periods P ~ 10~^""^ s and the mag- 
netic field strength B ~ iQio-^i^ G, and the old miUisec- 
ond radio pulsars, which have magnetic field strength 
as low as 10^"^^ G. While most radio pulsars are iso- 
lated objects, the millisecond pulsars are predominantly 
in binaries, suggesting that the magnetic fields decay 
with time, perhaps by an accretion of matte r from the 
binary companion fe.g.. IHarding fc Laill2006l and refer- 
ences therein). 

For the isolated radio pulsars, it is still an open ques- 
tion whether the NS magnetic fiel ds do or do not de - 
cay with time. For example, iNaravan fc Ostrikeii (|199O0 
argued that the field should decay exponentially on a 
few Myr timescale (see also ICordes fc Chernoffl 119981 : 
IGonthier et al.|[2004D . while s imilar studies implied the 
decay tim e > 100 Myr fe.g.. 'B hattacharva et al.l 119921 : 
iFauchcr-G iguere fc Kaspi 2006). Such divergent results 
may come from the uncertainties inherent to those pop- 
ulation synthesis studies, such as the selection effects, 
the luminosity evolution, and the dependence of beam- 
ing fraction on period. On the other hand, the discov- 
ery of magnetars has provided several important evi- 
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deuces to favor the i nagnetic field decay on relatively 
short timescales fe.g.. lArras et al.l 120041 : iHarding fc Lail 
|2006[) . In fact, the soft gamma-ray repeaters and anoma- 
lous X-ray pulsars, observed as young (< 10^ yr) NSs 
with strong (~ iQi^^io q-j magnetic fields, are consid- 
ered to be powered by the decay of their magnetic fields 
(e.g.. iWoods fc ThompsonI |2006[) . Giant flaring activi- 
ties are proposed to be the rapid r elease of the magnetic 
stress build ing in the NS crust () Thompson fc Duncan 
0995, 199^ ILvutikovl [20061) . More recently, |Eonset_aL 
pOOTf ) presented an evidence that the magnetic field de- 
cay of ~ 10 Myr, can explain the thermal evolution from 
magnetars continuously to ordinary radio pulsars. At 
present, these ideas and new observations seem to favor 
the existence of the magnetic field decay in the isolated 
NSs. 

The pioneering study by IGoldreich fc Reiseneggeii 
(|l992[ ) identified the dissipation processes of the mag- 
netic energy in the crust of the isolated NS during its 
evolution. They first showed that magnetized turbulence 
in the crust of NSs can be described by the electron mag- 
netohydrodynamic (EMHD) equation, in which the time 
evolution of the magnetic field is governed by the advec- 
tion of the field by the Hall drift, the Ohmic dissipation, 
and the ambipolar diffusion. On top of the Ohmic de- 
cay and the ambipolar diffusion, they first proposed that 
the Hall drift, though non-dissipative itself, could be an 
important ingredient for the field decay because it can 
lead to dissipation through the whistler cascade of the 
turbulence. As the eddy size of the turbulence becomes 
smaller, the magnetic energy was transported from large 
to small scales, leading to the dissipation of the magnetic 
energy finally via the Ohmic dissipation. 

In order to confirm their prediction, numerical simu- 
lations in EMHD are required because such a cascading 
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process is essentially a non-linear process. In the two- 
and three-dimensi onal (2D and 3D ) simulations of the 
low /3 plasma, Bisk amp et al.l ()1999[ ) showed a clear cas- 
cade of the energy due to the whistler turbulence over 
more than an order of mag nitude in the leng t h scale. The 
3D EMHD simulations bv lCho fc LazarianI (|2004l . |2009( ) 
confirmed the scale-dependent anisotropy in the EMHD 
turbulence. They also showed that the anisotropic cas- 
cading processes in the EMHD proceed via the propa- 
gation of the whistler waves along t he global magnetic 
fields (see, also lNarita fc Garvll201Cl[ ). This nature is in 
contrast to the conventional MHD turbulence, in which 
the Alfven waves play a majo r role to transport the tur- 
bulent magnetic energy ( e .g.. iGoldreich fc Sridharl 119951 : 
IBiskamp fc Wehe^ [l98l ICho et al.l \200W . Reflecting 
the fact that the propagation of the whistler waves are 
more dispersive than the Alfven waves, the power spec- 
tral density in the EMHD turbulence was shown to be- 
come steeper compared to the MHD turbulence. Here it 
should be noted in the EMHD simulations that one can 
precisely follow the cascading process, however the dis- 
sipation process should be treated phenomenologically, 
i.e., via the resistivity or th e so-called hyperdiffusivity 
(|Cho fc Lazarianl[200l [20091) . 

In this paper, we perform the Particle-In-Cell (PIC) 
simulations aiming to understand the decay process of 
the magnetic fields in NSs. As well known, the PIC 
simulations, which have been often p erformed to study 
the turbulent cascade in solar wind ()Saito et al.l 120081 : 
IGarv et al.l[2008l[20Tol) . can precisely take into account 
the electromagnetic modes in plasma by solving the full 
Maxwell equations coupling with the all species of plasma 
particles. The ordinary PIC simulation can treat the 
coUisionless plasma in which the coUisional frequency 
is smaller than the typical frequencies of plasma par- 
ticles, s uch as the plasma frequenc y and the gyro fre- 
quency. IShternin fc YakovlevI (j2006D have estimated the 
coUisional frequency between the particles in NSs. From 
their results, the electron coUisional frequency in the 
outer crust of the NS is of the order of 10~^ujpe with 
ujpe being the electron plasma frequency, which indicates 
that the electrons are marginally coUisionless in the elec- 
tron inertial scale. Since the gyro frequency is of the 
order of the plasma frequency, we approximately treat 
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Fig. 2. — Time evolution of the energy deviation. Blue, red, black 
curves are respectively for e = 0.8,0.4,0.1 and the solid, dashed, 
dot-dashed curves denote the electron kinetic energy (<5-Ekin)i the 
minus of the turbulent magnetic field energy whose field component 
is perpendicular to Bo {—SEtB,xz), and the minus of the turbu- 
lent magnetic field energy whose field component is parallel to Bq 
{-5EtB,y), respectively. 



the plasma as coUisionless to study the decay process of 
the magnetic fields in the outer crust of the NS. As a 
complement to the f oregoing EMHD simul ations (e.g., 
IBiskamp etall 119991 : iCho fc Lazaria^ 12004 [2009 ) that 
focus on the whistler cascade of the turbulence, the PIC 
simulations that we will present in this paper, should 
merit for understanding the mechanism of the magnetic 
field decay, cascading further down to the electron iner- 
tial scale. 

This paper is organized as follows. In Section 2, we 
show the numerical setup of the PIC simulations. The 
numerical results of the turbulent cascade and the dis- 
sipation process are shown in § 3. We summarize our 
results and discuss their implications in Section 4. 

2. NUMERICAL METHODS AND INITIAL CONDITIONS 

We carry out 2-dimensional fully rela tivistic PIC simu- 
lations (e.g.. lBirdsall fc Langdonll200lD . The basic equa- 
tions are described as. 
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Fig. 3. — (Upper panel) Power spectral density of the turbulent magnetic fields, log|Bt/iJoP at tuip^ = 1350 for e = 0.1,0.4,0.8 from 
left to right, respectively. The lower panels show the 1-dimensional, cumulative power spectral density. Thick and thin solid curves show 
\Bt{kx)/ Bq\^ and \Bt{ky)/ Bq]'^ , respectively, while dashed curves do the initial spectra of the turbulent magnetic energy. 
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where Pj, qi, E, B, pe, and j are the momentum, the 
particle charge, the electric fields, the magnetic fields, 
the charge density, and the current density, respectively. 
/3j = Vi jc is the three velocity normalized by the speed of 
light. The subscript i denotes the particle species. Here 
we ignore the quantum effects such as the Landau level, 
for simplicity. The particle motion is determined by solv- 
ing the special relativistic equation of motion. Then the 
electric charge and the charge current are obtained from 
equations ^ and ([7]), in which the shape factor S can 
extrapolate the physical quantities on discretized grids 
at X from particles. By using the electric charge and 
the current density, the electromagnetic fields are up- 
dated by solving Maxwell equations, so that the system 
can be solved self-consistently. We solve these equations 
in the rectangular coordinates by assuming dj [pz) = 



(the so-called 2.5-dimension). Number of grid points is 
{Nx,Ny) — (1024,1024) and the corresponding system 

size is Lx — Ly — 102 Ac/ojpe, where ujpe = Aime^ / rrie 
is the electron plasma frequency. Each cell contains 80 
particles and the total number of particles in the simula- 
tion box is 8.4 X 10^ for each species. The boundary con- 
ditions are periodic in both directions. We assume that 
ions are immovable, which is reasonable in the outer crust 
of the NS (e.g., Goldrcich & Rcisenegger 1992). We con- 
firmed that the numerical results are qualitatively and 
quantitatively consistent with those when we take into 
account the ion motion. 

In the initial state, the plasma distributes uniformly 
in space and the electron distribution function obeys the 
relat ivistic Maxwelli an with the temperature Tg = 10® K 
(e.g. JAguilera et al.>200 8'). We assume that there are two 
components of the magnetic fields in the NS. The first 
one is the uniform magnetic fields that would globally 
penetrate the NS, which should be necessary to explain 
the spin-down of the NS. The global field is set to be 
uniform in space B — BoBy, where By is the unit vector 
in y-direction. The corresponding ratio of the gyro fre- 
quency evaluated from the uniform magnetic fields ujce = 
eBo/{mec) and the plasma frequency ujpe is 0.5. The 
second one is the turbulent magnetic fields, which would 
be cascaded down from large to small scales due to the 
Hall turbulence. The initial turbulent magnetic fields are 
modeled to take the form of B{k) oc exp[—k^ /kg + iuk], 
where k is the wave number, ko = O.Slwpe/c is the typ- 
ical wave number of the initial fluctuations and ak is 
the random phase. The hat denotes the variable in the 
Fourier space. We change the amplitude of the turbulent 
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Fig. 4. — Time evolution of the anisotropy tan^ 9 for e = 0.8 
(Thick sohd), e = 0.4 (thin solid), and e = 0.1 (dotted curve). 



magnetic fields Bt — B B^ey by introducing a param- 
eter, e = / dF Bt{t = 0)7 / dV = 0.1,0.4,0.8. We 
also carry out the simulation with e = 0.0 to assess the 
validity of our simulation codes. 

3. RESULTS 
3.1. Energy Exchange 

Figure [T] shows the turbulent magnetic field energy 
density (color) and the magnetic field lines (white curves) 
at the initial state ttOpe = (left) and at the end of the 
simulation tujpe = 1350 (right) for e = 0.1. It is clearly 
shown that the turbulent magnetic field energy given at 
the initial state (bright spots in the left panel) decreases 
at the end of the simulations (right panel). The pattern 
size of the turbulent magnetic fields becomes relatively 
smaller with time. These results imply that the magnetic 
energy is converted to the plasma energy. 

To see clearly how the energy conversion proceeds, 
we calculate deviations of each energy from the ini- 
tial state, such as the electron kinetic energy ^E'kin = 
^jm[7(t) — j{t = 0)]/Eo and the turbulent magnetic 
field energies 6EtB,.. = J dV {{Bl + Bj) - [Bl{t = 
0)+Bl{t = Q)\}/{%ttEq) &nd5EtB,y = / dF [B^-Blit = 
0)]/{8ttEo), where Eo = f dVBf/{8TT) is the magnetic 
energy of uniform fields. Figure |2] depicts the time evo- 
lution of the energy deviation of the electrons JE'kin (solid 
curves) and of the minus of the turbulent magnetic ener- 
gies —SEtB,xz (dashed curves) and —SEtB,y (dot-dashed 
curves). It can be seen that the turbulent magnetic en- 
ergy in the perpendicular direction {SEtB,xz) decreases 
with time (note the minus sign), while the particle ki- 
netic energy ((5i?kin) increases with it. The energy de- 
viation of the particles (5i?kin is almost comparable to 
—5EtB,xz- This means that the perpendicular compo- 
nent of the magnetic energy is preferentially converted 
to the particle energy. Here it should be noted that the 
small difference between (5i?kin and —SEb.xz comes from 
the electric field energy generated by turbulent motions. 
The particle energy distribution function obtained in this 
simulation can be well fitted by the Maxwellian distribu- 
tion. Therefore, the plasma is not accelerated, but is 
heated up through the particle- wave interaction. The re- 
sulting heating is expected to become larger for models 
with larger initial turbulent magnetic energy (see models 
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Fig. 5. — Time evolution of the electron temperature. Blue, red, 
black, and green curves are, respectively, for e = 0.8,0.4,0.1,0.0 
and solid and dashed curves are for Tyy and T^x, respectively. 
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Fig. 6. — Time evolution of the ratio of the electron temperature 
Tyy/Txx- Blue, red, black and green curves are e = 0.8, 0.4, 0.1, 0.0, 
respectively. 



with different e in Figure 2). The efficient energy conver- 
sion to the perpendicular component also suggests that 
the energy dissipation will proceed anisotropically, which 
acts to make the configuration of the magnetic-field lines 
approach to the uniform one, i.e., the potential field. 

3.2. Anisotropy 

To see the anisotropy more in detail, we perform the 
Fourier analysis of the magnetic fields. Top panels of 
Figure [3] show the power spectral density (PSD) of the 

turbulent mag netic field |Bt/BoP at iwpe = 1350. It 
can be seen that the anisotropic turbulence preferen- 
tially transports the magnetic energy perpendicular to 
the uniform magnetic fields. The PSD is larger for a 
larger e, suggesting that the larger amplitudes of the tur- 
bulence lead to an efficient electron heating. The lower 
panels of Figure |3] show 1-dimensional cumulative power 
spectrum densities of the turbulent magnetic fields at 
tujpf, = 1350. The cumulative spectrum is defined as 

\Bt{kx)\'^ = Xlfc \Bt{kx,ky)\'^ (thick solid curves) and 
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Fig. 7. — Deviation of the electron distribution function at tuiri 



1350 from the initial state at tuirj 



denote the direction of the uniform magnetic fields Bq. Left and right panels denote for e = 0.1 and O.i 



= in /3a; — /3y plane. Blue arrows 
respectively. 



\Bt{ky)\'^ = J2k \Bt{kx:ky)\'^ (thin sohd). The summa- 
tion in fc-space is performed in the rage of < fc^^, fcy < 10 
to reduce thermal noises. Dashed curves show the initial 
turbulent spectrum, which is isotropic in k^-ky space. 
From these panels, the turbulent magnetic fields are 
shown to proceed via the forward cascade perpendicu- 
lar to the uniform magnetic fields. The PSD in the per- 
pendicular direction has a clear power-law distribution. 
Also the PSD in parallel direction also deviates from the 
initial (gaussian) distribution for a larger k. It suggests 
that the turbulent energy is mainly transported in per- 
pendicular direction, but a some part of the energy (a 
few percent) is transported in the parallel direction. 

To visualize the time evolution of the anisotropy, we 
show the time evolution of a quantity 6 in Figure |4l 



tan^ e 
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(see IShebalin et al.lll983HSaito et al.ll2008[ ). As already 
mentioned, the turbulent eddies interact each other with 
time more frequently for a larger e, transferring the tur- 
bulent energy to smaller scales. When tujpe ^ 100, the 
anisotropic cascades proceed faster for a larger e because 
the cascade rate is an increasing function of the fluctu- 
ation. The turbulent energy is then transferred to the 
smaller scale. After that, the anisotropy 9 saturates al- 
most at a constant level. The saturation level is shown 
to be highest for e = 0.4 among the computed models. 
This is because too much initial turbulent energy (like 
e = 0.8) disturbs the direction of the uniform magnetic 
fields, acting to smearing out the anisotropy. 

3.3. Electron Heating 

Now we are in a position to evaluate the tempera- 
ture of electrons heated by the particle- wave interactions 
mentioned above. In doing so, we utilize the energy- 
momentum tensor as 



me J d^p f{p) 



where e, and / are the electron energy, the electron 
four momentum, and the electron distribution function, 
respectively. It should be noted that the temperature 
obtained from this equation is evaluated in the observer 
frame, while it should be naturally defined in the comov- 
ing frame. However this simplification is good enough 
in our case because the averaged velocity of electrons in 
mostly random motions is much smaller than the speed 
of light. 

Figure [5] shows the time evolution of the electron tem- 
perature. It can be seen that the electron temperature 
parallel to the uniform magnetic field {Tyy, solid curves) 
increases with time faster than the perpendicular one 
{Txx^ dashed curves), and then Tyy becomes larger than 
Txx- It is also shown that the electrons are heated more 
rapidly for a larger e because the turbulent energy cas- 
cade proceeds faster. Figure [6] shows the time evolu- 
tion of the ratio of Tyy and Txx- Before tujpe < 200, 
the anisotropy of the electron temperature increases with 
time faster for a larger e. In this epoch, the electron is 
preferentially heated in the direction parallel to the uni- 
form magnetic fields. Later on, the ratio approaches to 
a nearly constant value. The saturation level is highest 
for e = 0.4 by the same reason as mentioned above (e.g.. 
Figure 111). 

Figure [7] shows the deviation of the electron distribu- 



tion function at tu)^ 



1350 from the initial state at 



(9) 



tujpe = for e = 0.1 (left panel) and e = 0.8 (right panel) 
on the plane of jSx — Py with /?; = Vi/c being the three 
velocity normalized by the speed of light. It can be seen 
that the distribution function is strongly distorted in the 
y-direction for e = 0.1. Considering the cascading pro- 
cess in perpendicular direction and the electron heating 
along the uniform magnetic fields, a Landau resonance 
{uj = kyVy) is the main mechanism for the electron heat- 
ing. The electromagnetic waves generated by the turbu- 
lent motion can interact with electrons, which results in 
the Landau resonance. 

Figure [5] shows the electron distribution function at 
tujpe = 1350 for e = 0.1 (left) and for e = 0.8 (right) on 
the plane of /3x — Py As mentioned above, the distri- 
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Fig. 8. — Electron distribution function at tujpe = 1350 for e = 0.1 



bution function is distorted in parallel direction because 
of the electron heating due to the Landau resonance. It 
should be noted that in both cases, the electrons with 
higher energies (i.e., a larger 13) are shown to get energies 
also in the direction perpendicular to the uniform mag- 
netic fields, thus making the distribution function almost 
isotropic (see also Figure [5] that also shows the decrease 
in the parallel component of the magnetic field). This 
indicates that not only the Landau resonance, but also 
the cyclotron resonance (a; = kyVy + oJce) plays a role 
of the electron heating and t he dissipation of the tur- 
bulent magnetic energy (e.g.. iSaito et al.|[2008D . While 
the direction of the electron heating is parallel to i?o for 
the Landau resonance, there is an energy exchange be- 
tween the waves and the perpendicular motion of the 
particles w hen the cyclotron resonance occurs (see, e.g., 
lGarvl[l993l ). It is also interesting to note that the en- 
ergy dissipation rate of the cyclotron resonance obtained 
in this simulation is almost comparable to that of the 
Landau resonance, which has been known to be the case 
only in the linear regime aMivamotoll2005D . The disper- 
sion relation of parallel propagating waves of Ex is shown 
in Figure ini (color contour). The dashed curves denote 
the dispersion relations for the i?-mode (top and bottom 
curves) and L-mode (middle curve) obtained from the 
linear analysis. Although the L-mode waves are gener- 
ated by the turbulent motion, it is noted that they do not 
exchange energies with electrons since their polarization 
is in an opposite direction to the mean magnetic fields. 
The L-mode electromagnetic waves can interact mainly 
with ions but they are now assumed to be immovable as 
already mentioned. Thus the turbulence is maintained 
by the i?-mode (whistler) waves in the electron scale. 

Comparing to the results of the linear analysis (dashed 
lines) , the dispersion relation obtained from the numeri- 
cal results is not so clear, mainly because the large ampli- 
tudes of the turbulent magnetic fields injected at the ini- 
tial state proceed forward cascades. By boldly extrapo- 
lating the agreement with the linear analysis with respect 
to the i?-mode wave (whistler wave) for kyc/wpe ^ 1.0 
(bottom line), the i?-mode wave in the higher ky regime 
would be damped due to the cyclotron resonance. The 
findings of this paper (i.e., the occurrence of the heating 
both in the parallel and perpendicular direction as well 




(left) and for e = 0.8 (right) on the plane of — /3i 



as the presence of the _R-mode damping) suggest that the 
Landau resonance and the cyclotron resonance would be 
pivotal factors for the dissipation of the magnetic energy 
and the resulting electron heating. 

4. SUMMARY AND DISCUSSIONS 

Bearing in mind the application to the outer crust of 
the NS, we investigated the turbulent magnetic field de- 
cay by means of the fully relativistic 2-dimensional PIC 
simulations. In the numerical simulations, the initial 
magnetic fields were set to be composed both of the uni- 
form magnetic fields that model the global fields pen- 
etrating the NS and of the turbulent magnetic fields 
that would be originated from the Hall cascade of the 
large-scale turbulence. Wc showed that the turbulent 
whistler cascade transports the magnetic energy pref- 
erentially in the direction perpendicular to the uniform 
magnetic fields. It was also found that the distribution 
function of electrons becomes anisotropic because elec- 
trons with lower energies are predominantly heated in 
the direction parallel to the uniform magnetic fields due 
to the Landau resonance, while electrons with higher 
energies are heated mainly by the cyclotron resonance 
that makes the distribution function isotropic for the 
high energy tails. Furthermore we pointed out that the 
degree of anisotropy takes maximum as a function of 
the initial turbulent magnetic energy. This is because 
too much initial turbulent energy disturbs the direction 
of the uniform magnetic fields, acting to smearing out 
the anisotropy. The findings of this paper suggest that 
the particle-wave interactions via the Landau resonance 
and the cyclotron resonance, alternative to the conven- 
tional particle-particle collisions via the Ohmic dissipa- 
tion, would be a pivotal dissipation mechanism in the 
outer crust of NSs. 

The initial turbulent magnetic fields assumed in this 
study are relatively larger conipared to the previous work 
(|Saito et al.ll200l iGarv et al.ll2008l ). For e < 0.4 in our 
simulation, the obtained results of the anisotropic cas- 
cade and the electron heating are basically consistent 
with the previous studies. As already mentioned, the 
anisotropy of the temperature decreases as e increases 
furthermore since the turbulent fields disturb the direc- 
tion of the uniform magnetic fields. Although the tem- 
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Fig. 9. — Dispersion relation of parallel propagating waves of 
Ex- The color shows the numerical result for e = 0.1. The dashed 
curves represent the dispersion relations for the /?-mode (top and 
bottom) and L-mode (middle) obtained from the linear analysis. 



perature becomes nearly isotropic, it should be noted 
that the cascading process itself is still anisotropic. Ac- 
cording to Cho & Lazarian (2004, 2009), the small eddy 
(with size ~ I) of the turbulent magnetic fields bi in- 
teracts not with the global magnetic fields Bq, but with 
the 'local' mean magnetic fields B^. Thus even when the 
turbulent magnetic field energy is much smaller than that 
of the uniform magnetic fields (especially Bq = 0), the 
cascading process should be locally anisotropic along B^. 
Such diffusion processes lead to the anisotropic heating 
of electrons, however the resulting electron temperature 
is in average isotropic, which should be the case obtained 
for high e in our simulation. 

In the conventional model of the magnetic field de- 
cay, the magnetic fields are considered to be dissipated 
through the collisional process, namely via the Ohmic 
dissipation. While in this paper, we proposed that the 
magnetic fields are dissipated due to the Landau and cy- 
clotron resonances, leading to the plasma heating in the 
collisionless regime. Since the wave-particle interactions 
occur in a very small scale (~ gyro radius), the corre- 
sponding timescale (~ w^^) is much shorter than that 



of the Ohmic dissipation, and it is instantaneous com- 
pared to the expected magnetic field decay timescale in 
the NSs. This means that the decay timescale should be 
determined by the cascading processes of the magnetic 
fields above the electron inertia scale. Although the PIC 
simulations have an advantage in its capability to de- 
termine the dissipation processes consistently, it is still 
computational too expensive to perform the PIC simu- 
lations covering over such a wide spacial range required 
to estimate the timescale. For the purpose, we think it 
important to perform the EMHD simulations including a 
phenomenological diffusivity which is adjusted to mimic 
the dissipation obtained in this study. Although this is 
apparently beyond the scope of this paper, we regard it 
as one of the most important tasks, which we plan to 
investigate as a sequel of this study. 

So far there have been extensive work that focuses on 
the origin of the o bserved anisotropy in the NSs' surface 
temperature (e.g.. iGeppert et al.ll2004 : iPons fc Gepperti 



[2557t lAguilera et al.ri2008l: iPons et al]|2009D . Unfortu 
nately it is hard for us to do so immediately because the 
anisotropy obtained in the current simulation is confined 
in a very small scale (collisionless scale). To clarify how 
the particle-particle collisions in the more larger scales 
redistribute the electron distribution function and the re- 
sulting global temperature could be, one may need a new 
simulation technique which bridges the PIC simulation 
and the global (E)MHD simulations, albeit not an easy 
job at present. In addition to the NSs' surface, we spec- 
ulate that the anisotropy due to the Landau resonance 
could be important also in the pulsar atmosphere because 
the plasma density is so small that the plasma there is 
expected to be globally collisionless. The consideration 
of the electron heating due to the Landau and cyclotron 
resonances might make the origin of the observed tem- 
perature anisotropy o f NSs (|Zavlinl [2007L iHaberll 120071 : 
iNakagawa et al.ll2009f) less mysterious, which we are go- 
ing to study one by one as an extension of this study. 



Numerical computations were carried out on Cray XT4 
at Center for Computational Astrophysics, CfCA, of Na- 
tional Astronomical Observatory of Japan and on Fujitsu 
FX-1 at JAXA Supercomputer System (JSS) of Japan 
Aerospace Exploration Agency (JAXA). This study was 
supported in part by the Grants-in-Aid for the Scientific 
Research from the Ministry of Education, Science and 
Culture of Japan (Nos. 19540309 and 20740150). 
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